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ABSTRACT 

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab initio quantum chem- 
istry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the 
ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic 
excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Mpller-Plesset perturbation 
theory (MP2), Z-averaged perturbation theory (ZAPT2), coupled cluster singles and doubles with perturbative triples corrections 
(CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ 
and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moder- 
ate (7-20 kcal moL 1 ) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between 
molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the ni- 
trogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet- 
visible spectra for electronic transitions with large oscillator strengths at the B3LYP, o>B97, and equations of motion coupled cluster 
(EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength elec- 
tronic transition in the visible region of the electromagnetic spectmm. 
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1. Introduction 

In the interstellar medium (ISM) thus far, inorganic salts, neutral 
hydrocarbons, olefins, cumulenes, nitriles, alcohols, carboxylic 
acids, cations, and anions have been identified (Gerakins 2011). 
Some of these are in the gas phase, and some are in interstellar 
ices and grains. Large polycyclic aromatic hydrocarbons (PAH) 
are also known to exist in the ISM, as a general class of interstel- 
lar molecules (Chillier et al. 1999; Salama et al. 1999). Recently, 
two fullerene molecules (C60 and C70) have been identified in a 
carbon rich planetary nebulae (Cami et al. 2010). It is proposed 
that the identification of individual PAH molecules could help 
to resolve some of the mysteries associated with the diffuse in- 
terstellar bands (DIBs; Salama et al. 1999; Snow 2001) and the 
galactic extended red emission (ERE; Rhee et al. 2007). Similar 
to PAHs, spectroscopic signatures of nitrogen containing poly- 
cyclic aromatic hydrocarbons (PANH) are prevalent in various 
cosmological environments (Hudgins et al. 2005). Nitrogen con- 
taining molecules and ions have been observed throughout the 
ISM. Gas phase reactions of nitrogen containing species with 
carbonaceous compounds represent an ideal way to incorporate 
nitrogen into carbon chain growth, and ultimately into cyclic and 
polycyclic carbon compounds. Although the smaller molecules 
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are pervasive in the ISM, how large molecules, chains, poly- 
cyclic complexes, and grains are formed in the ISM remains an 
open question. 

Weak intermolecular interactions dipole-dipole, ion- induced 
dipole, hydrogen bonding, and dispersion (van der Waals) forces 
between atoms, molecules, and ions, play major roles in nature. 
For example, in biological systems hydrogen bonding helps to 
keep a whole DNA strand together. In the gas phase, long-range 
ion-dipole and electrostatic interactions guide species to come 
together to form aggregates. In the very low particle density cold 
ISM, in the absence of other strong forces, these weak interac- 
tions can guide ions and molecules to condense into aggregates 
and initiate reactions (Herbst 2001; Smith 2002). Ions are rela- 
tively abundant species in the strong radiation fields of the ISM, 
in star forming regions, in planetary nebulae, and in boundary 
regions between hot and cold environments. Ion- molecule re- 
actions are barrierless and spontaneous in the gas phase. It is 
reasonable, therefore, to ask if this cold condensation can be a 
viable mechanism for forming small molecular complexes, large 
PAHs, and grains? And, if so, what are the specific pathways? 

Earlier theoretical and experimental studies have concen- 
trated on the n-n interactions leading to stable complexes in 
alkenes (Bates et al. 2008; Hopkins & Tschumper 2004), aro- 
matic systems (Sinnokrot & Sherrill 2004), and even in PAHs 
(Hirata et al. 1999; Szczepanski & Vala 1993; Weisman et al. 
2001; Weisman et al. 2005). The neutral complexes are typically 


Article published by EDP Sciences 


A74, page 1 of 12 


A&A 535, A74 (2011) 


bound by van der Waals interactions (though aromatic radical- 
radical stacking complexes with partial covalent character are 
known (e.g., Small et al. 2005). However, the cationic complexes 
are bound by both electrostatic and covalent interactions. Closed 
and open shell PAH dimers are estimated to be bound by ener- 
gies of about 20 kcal mol -1 , and absorb in the visible and some- 
times in the infrared (Rhee et al. 2007; Weisman et al. 2005). 
The cationic complexes of small olefins, especially the nitrogen 
containing ones, have not been explored in any detail. 

The objective of the current study is threefold. Firstly, we 
have examined, by means of ab initio quantum chemical meth- 
ods, the nature of binding among the simplest complexes be- 
tween some observed interstellar molecules and ions. In other 
words, the question of how the small precursors complex to- 
gether to form larger aggregates in interstellar low-density, low 
temperature conditions is explored. Secondly, we have computed 
the UV- Vis electronic spectra and the vibrational spectra of these 
species so that they may be compared with experiments as well 
as existing astronomical observations. It is vital to know the 
structural and spectroscopic nature of the ensuing complexes to 
understand their formation pathways. Some of these ion com- 
plexes may have large oscillator strengths in their electronic 
spectra due to charge transfer excitations. Thirdly, we have per- 
formed a systematic study of the complexes using various ab ini- 
tio quantum chemistry methods. Binding energies, excitation 
energies and geometrical parameters obtained using affordable 
density functional and perturbation theory methods were com- 
pared and benchmarked with results obtained by highly accurate 
coupled cluster methods. This benchmarking study will help us 
choose more affordable ab initio quantum chemistry methods for 
charting out pathways of formation of larger molecules in the 
future where coupled cluster calculations may not be computa- 
tionally feasible. 

Several studies have examined the radiative association re- 
action rates for relatively simple cation/neutral reactions such 
as C + + H 2 , and even up to systems with three heavy (non- 
hydrogen) atoms under ISM conditions (Black et al. 1975; Smith 
1992; Gerlich & Horning 1992; Barinovs & van Hemert 2006; 
Herbst et al. 2010). While our study is related to these investi- 
gations, examining formation products and their spectroscopic 
signatures, we have not investigated the radiative association 
rates, though we note that the studies cited above have demon- 
strated that in general radiative association rates increase as the 
cations/neutrals increase in size, which of course is consistent 
with the larger number of ro vibrational degrees of freedom. 

2. Theoretical methods 

The geometrical structures reported in this research are op- 
timized using second order Mpller-Plesset perturbation the- 
ory (MP2) along with Dunnings correlation consistent polar- 
ized valence triple zeta basis set (cc-pVTZ) and the augmented 
correlation consistent polarized valence triple zeta basis set 
(aug-cc-pVTZ; Dunning 1989). MP2/aug-cc-pVTZ optimized 
geometries were used when calculating the binding energies of 
the complexes employing the coupled cluster singles and dou- 
bles with a perturbative triples correction, CCSD(T), method 
(Purvis & Bartlett 1982; Raghavachari et al. 1989; Lee & 
Scuseria 1995), in conjunction with the cc-pVTZ and aug-cc- 
pVTZ basis sets. The binding energies were calculated by sub- 
tracting the absolute energies of the neutral and the cationic frag- 
ments from the energy of the cationic bound complex. Binding 
energies were also determined using a density functional the- 
ory (DFT) method for comparison. We employ Becke’ s widely 


used three-parameter hybrid functional (Becke 1993), based on 
the Becke exchange functional (Becke 1988) paired up with the 
correlation functional of Lee et al. (LYP; Lee et al. 1988). All 
geometries were fully optimized using B3LYP along with the 
cc-pVTZ basis set. A fine grid of 75 radial points and 590 an- 
gular functions was used for all the DFT computations. Open 
shell molecules were treated with unrestricted MP2 (UMP2), 
second order Z-averaged perturbation theory (ZAPT2) (Lee & 
Jayatilaka 1993; Lee et al. 1994), and UB3LYP methods along 
with the cc-pVTZ basis set. The ZAPT2 method is based on 
the use of a symmetric spin orbital basis (Jayatilaka & Lee 
1992), and leads to S x eigenfunctions in the n-particle space. 
It has been shown to have less spin contamination relative to 
UMP2 and other restricted open- shell perturbation theories, and 
to have a larger radius of convergence, and to be less costly 
(computationally) as well (Lee et al. 1994). For electronic exci- 
tation energies, starting from the ground electronic state to var- 
ious singlet excited states, including oscillator strengths, time 
dependent density functional theory (TDDFT) (Hirata et al. 
1999) computations were carried out using the B3LYP den- 
sity functional and the cc-pVTZ basis set. A range separated 
density functional, coB91 (Chai & Head-Gordon 2008), that 
is specifically designed for weak interactions was also em- 
ployed. For valence electronic transitions TDDFT has been 
shown to work well (Hirata et al. 1999; Weisman et al. 2001; 
Halasinski et al. 2003; Hirata et al. 2003). The range- separated 
functional, co B97, significantly reduces the serious problems of 
TDDFT with conventional functionals for charge-transfer ex- 
cited states (Dreuw & Head-Gordon 2004). The Q-Chem3.1 
(Shao et al. 2006) ab initio quantum chemistry code was used 
for all the DFT, TDDFT, MP2, (and UMP2) and ZAPT2 compu- 
tations. Coupled cluster calculations were carried out using the 
MOLPRO 2006 (Werner & Knowels 2006) quantum chemistry 
package. Harmonic vibrational frequency calculations were car- 
ried out at the MP2/cc-pVTZ and MP2/aug-cc-pVTZ levels to 
ensure that the molecules and complexes presented here repre- 
sent true minima. The Cartesian coordinates of all the species 
discussed in this work are included in the supplementary mate- 
rial available with the web site. 

3. Results 

Elementary growth processes starting from some of the nitrogen 
containing molecules or ions observed in interstellar space are 
investigated in the following sections. Specifically, several nitro- 
gen containing interstellar molecules and cations that have been 
observed in the ISM were paired up with small, unsaturated hy- 
drocarbons to simulate the initial nucleation step that could then 
be a pathway to formation of larger molecules, such as PAHs 
and PANHs, and ultimately to grain formation. Cationic com- 
plex formation is guided by long-range weak electrostatic ion- 
dipole type attraction, in which a wide variety of exotic products 
are expected. Some ion-molecule reactions in the low tempera- 
ture and ultra low-density conditions of the ISM readily make 
strong covalent bonds to produce linear or branched chains, 
replication of which leads to chain growth. In the next section 
we will see that these types of reactions take place with signif- 
icant favorable energy changes for larger molecule formation. 
Other dimers make relatively weak and chemically interesting 
face-to-face n - stacked or T-shaped complexes bound by electro- 
static or ion-dipole type interactions. Although several associa- 
tion or encounter complexes are investigated in depth for each 
system, only the simplest kinds of encounter complexes are in- 
cluded in this work as it is not intended to be an extensive search 
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Table 1 . Binding energies (in kcal mol ! ) of the dimer cations between (mostly) known interstellar molecules and cations are listed. 


Figure 
Fig. 1.A 
Fig. 1.B 
Fig. l.C 
Fig. l.D 
Fig. 2.A 
Fig. 2.B 
Fig. 3.A 
Fig. 3.B 
Fig. 3.C 
Fig. 4 
Fig. 5. A 
Fig. 5.B 
Fig. 6.A 
Fig. 6.B 


Reaction/Complex 
C 2 H 4 +HCNH + ■ 
C 2 H 4 +HCNH + ■ 
C 2 H 4 +HCNH + ■ 
C 2 H 4 +HCNH + ■ 

c 2 h 4 + h 2 cnhj - 
c 2 h 4 + h 2 cnhj - 

HCCH + HCNH + ■ 
HCCH + HCNH+ ■ 
HCCH + HCNH + ■ 


B 3LYP/cc-pVTZ MP2/cc-pVTZ MP2/aug-cc-pVTZ CCSD(T)/cc-pVTZ CCSD(T)/aug-cc-pVTZ 


[C 2 H4:HCNH] + 
[C 2 H4:HNCH] + 
[C 2 H 4 :HCNH] + , Boat 
[C 2 H 4 :HCNH] + , Bridge 
[C 2 H 4 :H 2 CNH 2 ] + 
[C 2 H 4 :H 2 CNH 2 ] + , Boat 
[HCCH:HCNH] + 
[HCCH:HCNH] + 
[HCCH:HCNH] + , Boat 
HCN+C 2 H+ -» [HCN:C 2 H 3 ] + 

HCC + + HCN — HCCHCN, Branched 
HCC + + NCH — [HCCNCH] + , Linear 
HCCC + + HCN -» [HCCC(H)CN] + , Branched 
HCCC + + HCN — [HCCCNCH] + , Linear 


Table 2. Binding energies (in kcal mol 


Figure Reaction/Complex 
Fig. 7 H 2 CCH 2 + H 2 CNH + - 
Fig. 8.A H 2 CNH 2 + H 2 CNHJ ■ 
Fig. 8.B H 2 CNH 2 + H 2 CNH 2 - 


> [H 2 CCH 2 :H 2 CNH] + 

> [H 2 CNH 2 ]+, Pi-stack 

> [H 2 CNH 2 ]+, Trans 


11.52 

11.97 

13.25 

11.25 

11.7 

16.34 

17.92 

19.71 

16.05 

16.58 

48.28 

50.88 

52.2 

49.47 

50.07 

7.88 

1.66 

3.16 

1.82 

1.96 

8.16 

8.3 

9.36 

6.89 

7.7 

26.00 

36.11 

38.85 

31.1 

32.01 

11.49 

8.51 

12.56 

11.06 

11.53 

15.93 

10.62 

18.49 

15.55 

16.08 

19.13 

15.7 

16.68 

20.28 

20.75 

72.01 

66.2 

71.57 

68.81 

68.38 

125.93 

125.62 

134.41 

154.13 

154.3 

180.73 

181.66 

198.26 

205.52 

205.6 

42.09 

38.25 

40.42 

38.48 

38.7 

59.92 

53.12 

55.03 

48.44 

48.81 

l complexes of astrochemically observed species. 



LYP/cc-pVTZ 

UMP2/cc-pVTZ 

UMP2/aug-cc-pVTZ 

ZAPT2/cc-pVTZ 

CCSD(T)/cc- 

17.31 

10.15 

11.25 

9.44 

10.28 

25.49 

22.84 

22.64 

23.71 

24.98 

35.76 

37.95 

38.42 

39.22 

40.2 


of the entire potential energy surfaces for all possible minima of 
the molecular complexes. That is, complex or molecule forma- 
tion in which a significant potential barrier must be overcome 
are not investigated in the present study. Binding energies for 
closed- shell complexes are presented in Table 1 and similarly 
for open-shell complexes in Table 2. In Tables 3 and 4, the elec- 
tronic excitation energies are presented at the B3LYP/cc-pVTZ 
and <uB97/cc-pVTZ levels of theories. In Table 3, the excitation 
energies of the complexes are given, while in Table 4 those of 
the monomeric molecules and ions are listed. While this paper 
is focused on nitrogen containing complexes, we have addition- 
ally studied the pure carbon complexes of HCCH, H 2 CCH 2 with 
their cations, which will be presented in another paper (Bera 
et al., in prep.). 

3.1. Closed shell cation 

I. [C2H 4 -HCNH] + . Ethylene, hydrogen cyanide and protonated 
hydrogen cyanides are observed species in some carbonaceous 
interstellar clouds (Snyder & Buhl 1971; Ziurys & Turner 1986). 
Electrostatic interaction between protonated hydrogen cyanide 
and ethylene leads to two C 2 V symmetry stable T-shaped com- 
plexes and two cyclic isomers with favorable binding ener- 
gies. The CH hydrogen in HCNH + interacts with the n elec- 
tron clouds of H 2 CCH 2 to make the stable complex A presented 
in Fig. l.A. The distance between the H atom and either of 
the C atoms of ethylene is 2.156, and the binding energy is 
13.24 kcal mol -1 at the geometry optimized at the MP2/aug-cc- 
pVTZ level of theory. As shown in Table 1 , the binding energy 
at the higher CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level 
of theory is 11.7 kcal mol -1 (henceforth we will abbreviate this 
as the CCSD(T)/aug-cc-pVTZ binding energy). The NH hydro- 
gen in HCNH + also interacts favorably with ethylene to make 
a more stable isomer, B, of this complex, presented in Fig. l.B. 
The NH-C distance is 1 .912 at the MP2/aug-cc-pVTZ optimized 
geometry. The more tightly bound nature of the complex is re- 
flected in the larger binding energy, 16.58 kcal mol -1 at the 
CCSD(T)/aug-cc-pVTZ level, as presented in Table 1. 


Although, formation of the lower energy isomer is expected 
due to thermodynamic reasons, the branching ratio for formation 
of the two isomers will depend on the approach geometries of 
the two constituents, and their relative translational energy. The 
apparent branching ratio may also depend upon the barrier to 
internal conversion between the isomers. 

The electronic excitation energies (in eV) and oscillator 
strengths are presented in Table 3 for all the complexes cal- 
culated using TDDFT as stated previously. Both have very 
strong oscillator strength absorption features in the UV. As seen 
from Table 3, B3LYP under predicts the excitation energies of 
these weakly bound complexes compared to the coB91 func- 
tional, which reduces self-interaction errors. The <x>B97 results 
(5.56 eV) are much closer to the excitation energies obtained at 
the EOM-CCSD/cc-pVTZ level of theory (6.46 eV) for the com- 
plex 1 A. For complex IB the EOM-CCSD excitation energy for 
the first excited state is 7.22 eV compared to 6.57 eV obtained 
with the tuB97 functional. 

Two other isomers with more traditional covalent bonding 
features are also located on the potential energy surface of the 
association complex. These are a four-membered ring and a 
three-membered ring structure. The four-membered ring isomer 
presented in Fig. l.C is bound by 50.07 kcal mol -1 energy at 
the CCSD(T)/aug-cc-pVTZ level, making it the most strongly 
bound of the isomers studied. By contrast, however, the three- 
membered ring structure (D) is the most weakly bound of the 
structures presented in Fig. l.D, with only 1.82 kcal mol -1 of 
binding energy. A comparison of the association energies for the 
complexes using the B3FYP, MP2 and CCSD(T) methods is pre- 
sented in Table 1. Both B3FYP and MP2 estimated association 
energies match very well in comparison to CCSD(T). The boat 
(Fig. l.C) and the bridge (Fig. l.D) isomers have intense elec- 
tronic absorptions at 7.99 eV and 8.54 eV respectively. 

Association energies for a complex between an ionic and a 
neutral species are generally large and thus are not going to be 
significantly changed by the basis set superposition error (BSSE) 
correction. One possible exception is the case of the bridge 
structure (Fig. l.D) for which the association energy is small 
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Table 3. TDDFT electronic excitation energies (in eV) of the complexes calculated with B3LYP/cc-pVTZ and o>B97/cc-pVTZ. 


B3LYP o>B97 


Figure 

Molecule 

El 

E2 

E3 

E4 

E5 

El 

E2 

E3 

E4 

E5 

Fig. 1.A 

C 2 H 4 + HCNH + -> [C 2 H4:HCNH] + 

3.81 

3.86 

5.64 

5.73 

7.14 

5.56 

5.73 

7.37 

7.98 

8.17 

l A ' 

Oscillator strength 

0.002 

0.00 

0.00 

0.00 

0.086 

0.014 

0.00 

0.00 

0.093 

0.199 

Fig. LB 

C 2 H 4 +HNCH + -» [C 2 H 4 :HNCH] + 

4.84 

4.91 

6.41 

6.59 

7.92 

6.57 

6.84 

7.62 

8.14 

8.27 

l A ' 

Oscillator strength 

0.004 

0.00 

0.00 

0.00 

0.00 

0.037 

0.00 

0.00 

0.00 

0.235 

Fig. l.C 

C 2 H 4 +HCNH + -> [C 2 H 4 :HCNH] + , Boat 

6.71 

7.41 

7.61 

7.93 

9.55 

7.01 

7.74 

7.99 

8.56 

9.79 

l A ' 

Oscillator strength 

0.00 

0.003 

0.105 

0.095 

0.001 

0.00 

0.002 

0.15 

0.076 

0.001 

Fig. l.D 

C 2 H 4 +HCNH + [C 2 H 4 :HCNH] + , Bridge 

3.51 

6.54 

6.7 

7.27 

7.97 

3.68 

6.98 

7.13 

7.65 

8.54 

l A' 

Oscillator strength 

0.001 

0.003 

0.029 

0.012 

0.248 

0.001 

0.041 

0.029 

0.025 

0.343 

Fig. 2.A 

C 2 H 4 +H 2 CNHJ [C 2 H 4 :H 2 CNH 2 ] + 

4.62 

5.86 

6.97 

7.37 

7.62 

5.32 

7.28 

7.61 

8.04 

8.44 

x A f 

Oscillator strength 

0.309 

0.00 

0.023 

0.139 

0.003 

0.314 

0.00 

0.245 

0.00 

0.001 

Fig. 2.B 

C 2 H 4 +// 2 CNHJ [C 2 H 4 :H 2 CNH 2 ] + , Boat 

8.81 

8.87 

9.66 

9.92 

10.03 

9.71 

9.97 

10.68 

10.87 

10.88 

X A 

Oscillator strength 

0.019 

0.004 

0.032 

0.004 

0.044 

0.026 

0.003 

0.031 

0.006 

0.005 

Fig. 3. A 

HCCH+HCNH + [HCCH:HCNH] + 

4.12 

4.31 

4.4 

4.5 

6.73 

5.8 

6.08 

6.39 

6.45 

7.12 

l A x 

Oscillator strength 

0.00 

0.016 

0.001 

0.00 

0.00 

0.00 

0.002 

0.00 

0.014 

0.00 

Fig. 3.B 

HCCH+HNCH + [HCCH:HNCH] + 

4.89 

5.1 

5.34 

5.45 

6.74 

6.19 

6.81 

7.39 

7.44 

7.51 

l A ! 

Oscillator strength 

0.00 

0.014 

0.001 

0.00 

0.00 

0.00 

0.001 

0.012 

0.00 

0.00 

Fig. 3.C 

HCCH+HNCH+ -> [HCCH:HNCH]\ Boat 

2.81 

5.62 

6.52 

7.37 

7.4 

3.03 

5.82 

6.85 

7.72 

7.92 

l A ' 

Oscillator strength 

0.005 

0.00 

0.00 

0.00 

0.027 

0.005 

0.00 

0.00 

0.00 

0.035 

Fig. 4 

HCN+C 2 H 3 + -> [HCN:C 2 H 3 ] + 

5.1 

5.76 

7.25 

7.65 

8.35 

5.73 

6.08 

7.52 

8.19 

8.68 

l A' 

Oscillator strength 

0.00 

0.3 

0.00 

0.00 

0.11 

0.00 

0.327 

0.00 

0.00 

0.264 

Fig. 5.A 

HCC + + HCN -> HCCHCN + , Branched 

3.06 

3.33 

4.9 

5.34 

6.39 

3.17 

3.39 

5.6 

6.04 

7.69 

l A' 

Oscillator strength 

0.00 

0.007 

0.011 

0.021 

0.365 

0.00 

0.009 

0.016 

0.177 

0.299 

Fig. 5.B 

HCC + + NCH [HCCNCH] + , Linear 

4.32 

4.46 

4.47 

7.46 

8.74 

4.46 

4.66 

4.67 

7.89 

8.67 


Oscillator strength 

0.00 

0.00 

0.00 

1.058 

0.00 

0.00 

0.00 

0.00 

1.249 

0.00 

Fig. 6.A 

HCCC + + HCN [HCCC(H)NCH] + , Branched 

1.62 

1.7 

4.32 

4.39 

5.11 

1.66 

2.26 

5.08 

6.41 

6.44 

l A ' 

Oscillator strength 

0.00 

0.003 

0.006 

0.001 

0.00 

0.00 

0.006 

0.187 

0.104 

0.00 

Fig. 6.B 

HCCC + + HCN [HCCCNCH] + , Linear 

0.83 

3.21 

4.32 

4.62 

5.37 

0.89 

3.44 

4.54 

5.47 

5.74 

l A ' 

Oscillator strength 

0.001 

0.001 

0.028 

0.002 

0.315 

0.001 

0.001 

0.062 

0.003 

0.343 

Fig. 7 

H 2 CCH 2 +H 2 CNH + [H 2 CCH 2 :H 2 CNH] + 

2.27 

3.39 

3.81 

4.15 

4.4 

2.29 

3.67 

4.09 

4.61 

4.77 

2 A ' 

Oscillator strength 

0.23 

0.00 

0.001 

0.003 

1.00 

0.21 

0.002 

0.003 

0.001 

0.01 

Fig. 8.A 

H 2 CNH 2 + H 2 CNH+ [H 2 CNH 2 ]+, Pi-stack 

2.86 

4.29 

4.94 

5.16 

5.23 

2.5 

4.89 

5.43 

5.48 

5.78 

2 A ' 

Oscillator strength 

0.208 

0.018 

0.00 

0.00 

0.00 

0.237 

0.008 

0.015 

0.00 

0.00 

Fig. 8.B 

H 2 CNH 2 + H 2 CNHJ [H 2 CNH 2 ]+, Trans 

3.01 

4.99 

5.68 

6.16 

6.38 

3.02 

4.99 

5.68 

6.16 

6.38 

2 A ' 

Oscillator strength 

0.195 

0.00 

0.00 

0.002 

0.00 

0.195 

0.00 

0.00 

0.002 

0.00 


Notes. Underneath the excitation energies the dimensionless oscillator strengths are presented. The electronic ground state is noted in the column 
with the figure name (the row below the figure name). 1 eV = 23.0605 kcal mol -1 = 1.16045 x 10 4 K. 


(1.96 kcal/mol). However, as Fig. l.D shows, after association, 
the geometry changes significantly from the initial monomer ge- 
ometries, which also explains why the association energy for an 
ion/neutral complex is small in this case in the association some 
new bonds are formed, but also some other bonds are weakened. 
A BSSE correction in the situation where the two monomer units 
do not retain their monomer quality but rather react to become 
one molecule is uninformative. Thus, BSSE corrections have not 
been computed for the ion/neutral complexes included in the 
present study. 


II. [H 2 CCH 2 -H 2 CNH 2 ] + . Methyleneimine, H 2 CNH, is ob- 
served in various interstellar environments (Johnson & Lovas 
1972). Planar protonated methyleneimine can bind with ethylene 
by a face-to-face stacking interaction. The parallel-displaced 
face stacked dimer cation of these two species is presented as 
isomer A in Fig. 2. The H^NHJ moiety is slipped from the cen- 
ter upon optimization (this complex has Cs symmetry), with the 
nitrogen pointing outwards at a right angle. At the MP2/aug-cc- 
pVTZ optimized geometry, the nearest carbon-carbon distance 
between the two moieties is 2.781, and the binding energy is 
9.36 kcal mol -1 relative to the ethylene and methylene amine 


cation dissociation channel. Upon further single-point refine- 
ment at the CCSD(T)/aug-cc-pVTZ level of theory, the bind- 
ing energy is 7.70 kcal mol -1 . This relatively weak dimer com- 
plex is bound by the attraction of the diffuse n electron cloud 
and the positive charge of the closed shell methylene amine 
cation. This example shows how simple double-bonded species, 
including nitrogen containing ones, can form stacked complexes 
with small but non-negligible binding energies by electrostatic 
interactions when one monomer is charged. The [H 2 CCH 2 - 
H 2 CNH 2 ] + complex, Fig. 2. A, has several singlet- singlet elec- 
tronic excitation bands in the UV; one around 4.62 eV (266 nm) 
has the strongest oscillator strength of 0.31. At the u>B97/cc- 
pVTZ level of theory the excitation energy for the same strong 
excitation band is 5.32 eV, which is much closer to the exci- 
tation energy, 5.62 eV, obtained at the EOM-CCSD/cc-pVTZ 
level of theory. The excitation energies calculated using B3LYP 
and qjB91 for all complexes studied are presented in Table 3. 
A more strongly bound four-membered cyclic structure, pre- 
sented as B in Fig. 2, is also located, with 31.10 kcal mol -1 
binding energy. This cyclic, slightly puckered structure is in- 
teresting as it incorporates a N atom in the ring, indicating the 
ease of ring formation in the gas phase due to electrostatic at- 
traction. Which of these two product channels is obtained will 
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Table 4. Excitation energies (in eV) and oscillator strengths of selected molecules and ions using time-dependent density functional theory 
(TDDFT) with B3LYP/cc-pVTZ and o>B97/cc-pVTZ. 


B3LYP o>B97 


Molecule 

El 

E2 

E3 

E4 

E5 

El 

E2 

E3 

E4 

E5 

HCCH, Excitation energy 

8.39 

8.67 

8.67 

12.42 

12.42 

8.34 

8.72 

8.72 

13.17 

13.17 

1 Z^, Oscillator strength 

0.00 

0.00 

0.00 

0.007 

0.007 

0.00 

0.00 

0.00 

0.009 

0.009 

H 2 CCH 2 , Excitation energy 

7.89 

7.91 

8.18 

8.6 

8.83 

7.89 

8.33 

8.6 

9.38 

9.64 

l A g , Oscillator strength 

0.363 

0.028 

0.00 

0.00 

0.00 

0.385 

0.00 

0.034 

0.00 

0.00 

HCN, Excitation energy 

7.71 

8.00 

8.00 

9.17 

9.17 

7.63 

8.02 

8.02 

9.41 

9.41 

1 E + , Oscillator strength 

0.00 

0.00 

0.00 

0.019 

0.019 

0.00 

0.00 

0.00 

0.019 

0.019 

HCNH + , Excitation energy 

8.39 

8.67 

8.67 

12.42 

12.43 

8.36 

8.71 

8.71 

13.17 

13.17 

1 E, Oscillator strength 

0.00 

0.00 

0.00 

0.007 

0.007 

0.00 

0.00 

0.00 

0.009 

0.009 

H 2 CNHJ, Excitation energy 

7.94 

9.52 

10.8 

11.33 

12.2 

8.04 

9.53 

10.91 

11.97 

12.79 

l A\, Oscillator strength 

0.00 

0.297 

0.008 

0.032 

0.00 

0.00 

0.302 

0.009 

0.039 

0.00 

H 2 CNH + , Excitation energy 

3.58 

5.59 

5.85 

6.89 

8.34 

3.49 

5.65 

5.91 

7.17 

8.58 

l A f , Oscillator strength 

0.00 

0.00 

0.001 

0.014 

0.006 

0.00 

0.00 

0.002 

0.013 

0.004 

H 2 CNH 2 , Excitation energy 

3.05 

4.22 

4.33 

6.13 

6.65 

3.89 

5.11 

5.11 

6.72 

7.13 

2 A\ Oscillator strength 

0.00 

0.00 

0.007 

0.00 

0.127 

0.00 

0.008 

0.00 

0.131 

0.00 


Notes. The electronic ground state is noted in the column with the molecule name. 



Fig. 1. MP2/aug-cc-pVTZ optimized geometry of the encounter com- 
plexes between HNCH + and H 2 CCH 2 . Bond lengths are expressed in A 
and angles in degrees. Binding energies (BE, in kcal mol -1 ) are pre- 
sented next to the figure labels. 


depend to some extent on the relative orientation of the reactant 
partners as they reach association distance, though there will be 
a large thermodynamic driving force towards isomer B . Further 
association of ring structure B with acetylene might produce a 
six-membered carbon/nitrogen cyclic molecule. 




Fig. 2. Electrostatically bound closed shell ^--complex and boat struc- 
ture between the H 2 CCH 2 and H 2 CNH 2 ion, presented here at the 
MP2/aug-cc-PVTZ level of theory. Bond angles are expressed in A and 
angles in degrees. Binding energies (BE, in kcal mol -1 ) are presented 
next to the figure labels. 


III. [HCCH-HCNH] + and [HCCH-HNCH] + . Acetylene and 
HCNH + are both observed in the ISM (Ziurys & Turner 
1986). In total, three isomers that can be formed by barrier- 
less association between HCCH and HCNH + have been located 
on the singlet potential energy surface of the complex. Two 
T-shaped weakly bound complexes as well as a slightly more 
strongly bound cyclic isomer are presented in Fig. 3. In both 
T-shaped structures, shown as A and B in Fig. 3, the terminal 
H atoms of HCNH + are attracted by the n electron cloud of 
HCCH. In HCNH + , the positive charge is distributed over all 
four atoms, but the N-terminal H atom has the most positive 
charge, hence the electrostatic attraction between that H atom 
and the n electron cloud of HCCH is stronger relative to the 
C-terminal H atom, and therefore the N-terminal H atom makes 
a stronger bond with the acetylene than the C terminal H-atom. 
The binding energies are presented in Table 1. The [HCCH- 
HCNH] + complex is bound by 12.56 kcal mol -1 at the MP2/aug- 
cc-pVTZ level and 11.53 kcal mol -1 at the CCSD(T)/aug-cc- 
pVTZ level with respect to the HCCH and HCNH + dissociation 
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1.0731 


1.1441 


1.0781 




Fig. 3. Electrostatically bound closed shell ^--complex and boat struc- 
ture between the H 2 CCH 2 and H^NHJ ion, presented here at the 
MP2/aug-cc-PVTZ level of theory. Bond angles are expressed in A and 
angles in degrees. Binding energies (BE, in kcal mol -1 ) are presented 
next to the figure labels. 


channel. Dissociating the more strongly bound [HCCH- 
HNCH] + T-shaped complex into HCCH and HCNH + will re- 
quire 18.49 and 16.08 kcal mol -1 at the MP2/aug-cc-pVTZ 
and CCSD(T)/aug-cc-pVTZ levels, respectively. As seen in 
Fig. 3, [HCCH-HNCH] + has a shorter NH distance (1.895) rel- 
ative to the CH distance, in [HCCH-HCNH] + , 2.110. Perhaps 
surprisingly, unlike the T-shaped complexes of HCNH + and 
H 2 CCH 2 discussed earlier (in Fig. 1), these complexes have 
weak electronic transitions in the ultraviolet. There is a cyclic 
isomer, presented in Fig. 3, with a four-membered ring and a 
slightly larger binding energy, relative to [HCCH-HNCH] + , of 
20.75 kcal mol -1 at the CCSD(T)/aug-cc-pVTZ level. This is 
the simplest cyclic structure that incorporates a N atom into 
its framework. Calculated electronic excitation energies for all 
structures are presented in the Table 3. The difference in exci- 
tation energies between B3LYP (4.12 eV) and a)B91 (5.80 eV) 
DFT methods stems from the fact that B3LYP is not suitable 
for predicting excitations that have an associated charge transfer 
component. Of the two density functionals employed here, the 
u>B97 excitation energy comes much closer to that obtained by 
the EOM-CCSD/cc-pVTZ method (7.0 eV). 

IV. [HCN-H 2 CCH] + . Another potentially viable mechanism for 
carbonaceous chain growth involves association of very abun- 
dant nitriles, such as HCN, and unsaturated molecular cations 
or fragment cations like vinyl cation. Steered by electrostatic 


Fig. 4. Electrostatically bound closed shell ^--complex and boat struc- 
ture between the H 2 CCH + and HCN, presented here at the MP2/aug- 
cc-PVTZ level of theory. Bond angles are expressed in A and angles 
in degrees. Binding energies (BE, in kcal mol -1 ) are presented under 
the figure. 


ion-dipole attraction, vinyl cation can be associated with polar 
hydrogen cyanide, incorporating nitrogen in the chain growth 
process. Association of HCN with vinyl cation results in the 
formation of a strong covalently bound closed shell molecu- 
lar cation. The MP2/aug-cc-pVTZ optimized geometry of this 
species, vinyl isocyanide cation, is shown in Fig. 4. Loss of a 
hydrogen from this complex leads to Cs^N 4- which was found 
to form in the selected ion flow tube during UV induced disso- 
ciation of pyrimidine and may be a key species in pyrimidine 
formation (Bera et al. 2009; Fondren et al. 2007). 

In vinyl isocyanide cation (see Fig. 4), the lone-pair elec- 
tron on nitrogen donates into the vacant carbon n orbital to 
form a strong CN sigma bond. The central sp 2 hybridized car- 
bon makes a short single bond, 1.388, with the sp hybridized 
nitrogen. At the MP2/aug-cc-pVTZ level of theory the bind- 
ing energy is 71.57 kcal mol -1 relative to the H 2 CQHL and 
HCN dissociation channel. At the CCSD(T)/aug-cc-pVTZ level, 
the binding energy is 68.38 kcal mol -1 using the geometries op- 
timized at the MP2/aug-cc-pVTZ level. The singlet ground state 
to singlet excited state vertical excitation energies obtained with 
TDDFT at the B3LYP/cc-pVTZ level of theory shows that the 
molecule has two relatively strong peaks in the UV at 8.35 eV 
(149 nm) and 5.76 eV (215 nm), with oscillator strengths of 0. 1 1 
and 0.30, respectively. The excitation energies (in eV) and oscil- 
lator strengths calculated at the B3LYP/cc-pVTZ and o>B97/cc- 
pVTZ levels of theories are included in Table 3, and are in quite 
good agreement with each other. 

V. [HCN-HCC] + . The neutral molecule acetylene (HCCH) and 
the acetyl radical (HCC) are observed species in the ISM, as 
is neutral C2. Protonated carbon CH + has also been observed, 
however, HCC + is as yet unobserved. Protonation of C2 will 
most likely produce HCC + in the ISM, though the abundance 


A74, page 6 of 12 


P. P. Bera et al.: Initiating molecular growth in the ISM 



A, BE 154.3 



Fig. 5. Closed shell complexes of HCC + and HCN or vice versa, pre- 
sented here at the MP2/aug-cc-PVTZ level of theory. Bond angles 
are expressed in A and angles in degrees. Binding energies (BE, in 
kcal mol -1 ) are presented under the figure. 

of HCC + in the ISM may be low due to its fast exothermic 
reaction with H 2 (McWan et al. 1999). Reaction between two 
closed shell species, acetyl cation (HCC + ) and hydrogen cyanide 
leads to two structural isomers: branched HCNC(H)C + and lin- 
ear HCCNCH + , whose MP2/aug-cc-pVTZ structures are pre- 
sented in Fig. 5. A n complex between closed shell HCC + 
and HCN is not stable due to very strong electrostatic attrac- 
tion. The N lone pair of electrons in HCN can donate into the 
empty n orbital of the acetyl cation to make a strong sigma 
bond. The binding energies for these two isomers are very 
strong, 134.42 and 198.26 kcal mol -1 , respectively at the same 
level of theory. Binding energies refined at the CCSD(T)/aug- 
cc-pVTZ level for the HCC + and HCN dissociation channel are 
154.30 and 205.59 kcal mol -1 for branched (Fig. 5) and linear 
isomers, respectively. The branched isomer has ten electronic 
absorption lines in the ultraviolet ranging from approximately 
8.59 eV to 3.06 eV (144 nm to 405 nm). One of these bands at 
6.39 eV (195 nm) exhibits a relatively strong oscillator strength 
of 0.36, while all others are weak. The linear isomer has only one 
very strong electronic excitation line at 166 nm with oscillator 
strength 1.1 computed at the TDDFT B3LYP/cc-pVTZ level of 
theory. Other excitations have nearly zero oscillator strength. For 
this case, the excitation energies and oscillator strengths com- 
puted with the B3LYP and range separated <x>B97 functionals 
are in good agreement with each other (Table 3). 

VI. [HCCC-HCN] + . The encounter between the HCCC + ion 
and HCN leads to two structural isomers, presented in Fig. 6. 
A n complex of HCCC + with HCN does not exist. However, the 
cr-bonded complexes are stable. The geometries optimized at the 
MP2/aug-cc-pVTZ level of theory are presented in Fig. 6. 
The binding energies, presented in Table 1, for the branched and 
the linear isomers are 40.42 and 55.03 kcal mol -1 respectively 
at the above mentioned level of theory. Binding energies relative 
to the HCCC + and HCN dissociation channel further refined at 
the CCSD(T)/aug-cc-pVTZ level are 38.70 and 48.8 1 kcal mol -1 



A, BE 38.7 



B, BE 48.81 

Fig. 6. Closed shell complexes of HCCC + and HCN or vice versa, 
presented here at the MP2/aug-cc-PVTZ level of theory. Bond angles 
are expressed in A and angles in degrees. Binding energies (BE, in 
kcal mol -1 ) are presented under the figure next to the labels. 


for the branched and the linear isomers (see Fig. 6). Both iso- 
mers have very strong electronic absorption signatures in the 
UV. One electronic excitation at 0.83 eV (B3LYP) or 0.89 eV 
(wB97), as presented in Table 3, falls into the infrared region. 
Oscillator strengths and excitation energies for these complexes 
computed using TDDFT B3LYP/cc-pVTZ and tuB97/cc-pVTZ 
are presented in Table 3. 

3.2. Open shell cations 

VII. [H 2 CCH 2 -H 2 CNH] + . The open shell cationic complex of 
H 2 CNH + (H 2 CNH is an observed species in the ISM) (Gerakins 
2011) with ethylene is an intriguing system. One proton linked 
to the carbon atom in H 2 CNH + is bound with the 7r-electron 
cloud of H 2 CCH 2 in this cationic complex. The binding en- 
ergy between H 2 CCH 2 and H 2 CNH + is 11.25 kcal mol -1 at 
the UMP2/aug-cc-pVTZ level (Table 2). The optimized geom- 
etry is presented in Fig. 7. ZAPT2/cc-pVTZ predicts the bind- 
ing energy to be 9.44 kcal/mol and at the CCSD(T)/cc-pVTZ 
level, the binding energy is 10.28 kcal mol -1 . The valence elec- 
tronic excitation energies and oscillator strengths computed with 
the B3LYP/cc-pVTZ and wB97/cc-pVTZ TDDFT methods are 
given in Table 3 and presented in Fig. 9 as a spectrum. The 
spectrum has multiple lines in the UV/VIS, only one of which 
at 2.27 eV (546 nm), has a relatively large oscillator strength 
of 0.23. Molecular orbital pictures indicate a charge transfer 
transition occurring from the highest doubly occupied molec- 
ular orbital (HDOMO) to the singly occupied molecular orbital 
(SOMO) (see Fig. 9) resulting in a very large oscillator strength 
at a very low wavelength for this electronic transition. For com- 
parison, the excitation energies of the monomeric species are 
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Fig. 7. Structure of weak electrostatically bound association complex 
of H 2 CNH + with H 2 CCH 2 molecule, presented here at the MP2/aug- 
cc-PVTZ level of theory. Bond angles are expressed in A and angles 
in degrees. Binding energies (BE, in kcal mol -1 ) are presented under 
the figure. 


also presented in Table 4. It is clear that the low energy, intense 
electronic excitation band is a feature that is exclusive to the 
complex as there is no such band for either C2H4 or H 2 CNH + . 

VIII. [H2CNH2 — H 2 CNH 2 ] + . Two isomers of the cationic dimer 
of H2CNH2 are presented in Fig. 8. One is a linear chain 
with carbon-carbon connectivity and the other is a n - stacked 
complex. The trans isomer, isomer A in Fig. 8, has an elon- 
gated 1.774 CC bond, and is bound by 38.42 kcal mol -1 at 
the UMP2/aug-cc-pVTZ level. Binding energies calculated us- 
ing ZAPT2 and CCSD(T) with cc-pVTZ basis set are 39.22 and 
40.20 kcal/mol respectively. The MP2/aug-cc-pVTZ optimized 
structure for isomer B (the 7r-stacked complex, see Fig. 8) shows 
that the carbon nitrogen distance between the monomer units 
is 2.622, which is an unusually long bond. The 7r-stacked com- 
plex is moderately strongly bound with a UMP2/aug-cc-pVTZ 
binding energy of 22.64 kcal mol -1 . The binding energies pre- 
dicted by the B3LYP (25.49) and ZAPT2 (23.7 1) levels of theory 
are also in agreement with that predicted by CCSD(T)/cc-pVTZ 
(24.98 kcal/mol). TDDFT calculations, Table 3, reveal electronic 



Fig. 8. The structures of the trans association complex and -stacked as- 
sociation complex of H 2 CNH 2 with H 2 CNHJ cation computed at the 
MP2/aug-cc-PVTZ level of theory. Bond angles are expressed in A and 
angles in degrees. Binding energies (BE, in kcal mol -1 ) are presented 
under the figure. 


excitation energies with significant oscillator strength for both 
isomers due to charge-transfer type transitions. The trans isomer 
has a strong excitation transition at 3.01 eV (411 nm) with an 
oscillator strength of 0.19. The 7r-stacked isomer has a strong 
electronic excitation transition with 0.21 oscillator strength at 
2.86 eV (433 nm). 

3.3. Excitation energies 

Excitation energies were calculated using time dependent den- 
sity functional theory with two different functionals, the widely 
used B3LYP and coB97 (Chai & Head-Gordon 2008), for com- 
parison purposes. The electronic excitation calculations reported 
here refer to bands that would be observed in absorption from 
the ground electronic state to an excited state. Emission bands, 
which would likely be similar but shifted, can be calculated, but 
for this work we have limited the scope to absorption bands. 
Excitation energies and oscillator strengths of the complexes 
are presented in Table 3, and that of the constituting molecules 
and ions are presented in Table 4. For tightly bound covalent 
complexes, the excitation energies predicted by both the density 
functionals are in good agreement with each other. However, 
for the non-covalently bound complexes the electronic excita- 
tion energies predicted by B3LYP and coB97 functionals can 
be different by more than 1 eV at times. As mentioned in the 
methods section, standard density functionals can perform very 
poorly for charge- transfer type excitations due to self-interaction 
errors, while range- separated hybrid density functionals, such as 
coB97, have been shown to perform well in predicting excita- 
tion energies of charge transfer states. We have performed more 
accurate equations of motion coupled cluster singles and dou- 
bles (EOM-CCSD) calculations for the cases where the B3LYP 
and co B97 excitation energies are found to be in disagreement. 
The EOM-CCSD calculations for the lowest five excited sates 
for the non-covalently bound complexes show very good agree- 
ment with those calculated using the coB97 density functional. 

We have computed the excitation energies of the five lowest 
energy excitations of the constituent neutral and ionic species 
and presented them in Table 4. A comparison of the excitation 
energies of the constituent monomer neutrals and ions with that 
of the complexes shows that none of the electronic excitations 
of the complexes are due to excitation in one constituent only. 
Some of the excitations are charge transfer in nature between the 
two moieties in the complex (e.g., Figs. l.A,B) while others are 
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due to complete reorganization of electronic structure after com- 
plexation. For example, first excitation energies of H 2 CCH 2 and 
H 2 CNH 2 , given in Table 4, are 7.89 eV and 8.04 eV respectively. 
The lowest excitation energy of the complex appears at much 
lower energy, 5.32 eV as given in Table 3, in the near ultra violet 
due to a charge transfer type transition between H 2 CCH 2 and 
H 2 CNH 2 . The electronic excitation energies produced by us- 
ing the hybrid functional B3LYP and the range separated u;B97 
agree well for all the species in Table 4 (unlike in the complexes 
where charge transfer plays a role). These excitation energies of 
the constituent monomer units can themselves be valuable for 
observation of these species in the ISM. 

3.4. Astrophysical implications 

The association complexes of mostly observed molecules and 
ions upon encounter are very interesting from an astrophysical 
perspective. As found in the present study (Tables 1 and 2), the 
ion-molecule complexes form with significant association en- 
ergies that must be dissipated. In the low-density of the ISM, 
collisions are rare, and therefore, this excess energy is not read- 
ily dissipated to third bodies. The excess energy is then either 
lost radiatively or retained in the internal degrees of freedom 
of the complexes making the complexes vibrationally and rota- 
tionally hot. The weakly bound complexes may therefore have 
very short lifetimes and dissociate or isomerize to give rise to 
species other than their parents. However, in the current study 
we have not investigated dissociation rates for different chan- 
nels. This could be investigated in the future. The open- shell 
cations have large oscillator strengths in the visible or near ul- 
traviolet region. Therefore, the open shell nitrogen containing 
complexes, such as the one presented in Fig. 7 with strong elec- 
tronic absorption lines in the visible (Fig. 9), are candidates for 
future laboratory study, and then observational detection. Most 
of the closed- shell association complexes have large oscillator 
strengths in the ultraviolet region of the electromagnetic spec- 
trum. Linear conjugated carbonaceous cumulenes are common 
species in the ISM. Therefore, the linear and branched closed 
shell association ion-complexes, such as the ones presented in 
Figs. 4-6, should also be good candidates for future observa- 
tion. From a growth perspective these nitrogen containing ionic 
association complexes can also serve as possible points for nu- 
cleation into larger molecules, such as nitrogeneous PAHs. The 
method of incorporation of nitrogen in the carbon framework is 
also related to how nitrogenous polycyclic molecules could form 
in the ISM. Identification of some of these intermediate com- 
plexes will help in understanding the involvement of nitrogen in 
the formation process of carbonaceous molecules in ISM. 

4. Conclusions 

Gas phase ion-dipole encounter complexes of observed, and 
some yet unobserved, interstellar molecules have been studied 
using density functional theory, Mpller-Plesset perturbation the- 
ory, and coupled cluster theory methods along with correla- 
tion consistent basis sets. Our main conclusion is that the sim- 
plest interstellar cations, steered by electrostatics, can initiate 
grain formation by forming larger chains, small cyclic cations, 
^-complexes, and strong ion-dipole complexes. In more detail, 
cations of some alkenes and alkynes, and nitrogen substituted 
alkenes and alkynes have been paired with their neutral coun- 
terparts making both closed shell and open shell complexes. 
Several different orientations, e.g. parallel, parallel displaced and 
T-shaped, of these encounter complexes have been tested for 


[H 2 CCH 2 -- H 2 CNHr'opensheH TDDFT/cc-pVT2 

- 
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Fig. 9. The absorption spectrum (oscillator strength vs. wave- 
length (nm)) of the H 2 CCH 2 :H 2 CNH + complex computed using time 
dependent density functional theory (TDDFT), at the B3LYP/cc- 
pVTZ level. The inset shows the highest doubly occupied molecular 
orbital (HDOMO) and the singly occupied molecular orbital (SOMO) 
of the complex that has a large oscillator strength. 


their stability, binding energy and electronic spectra. Acetylene 
and ethylene, when paired with protonated hydrogen cyanide 
(i.e. HCNH + ), all of which are abundant in the ISM, make strong 
T-shaped closed- shell encounter complexes in which the pro- 
ton plays a central role in bonding. Binding energies, ranging 
from 11.06 to 16.05 kcal mol -1 , in these complexes are too 
strong to be a conventional hydrogen bond and too weak to be 
a conventional cr-bond. These examples include the T-shaped 
complexes in Figs. 1 and 3. Ethylene and H 2 CNH 2 form a strong 
parallel-rotated n complex. Ethylene and H 2 CNH 2 also form a 
more stable four-membered cyclic structure, though its bind- 
ing energy is significantly less than a typical covalent bond. 
These examples are shown in Fig. 2. Closed shell n complexes 
of HCN and species with unfulfilled valencies (HCC + , HCCC + , 
and vinyl cation) are unstable upon geometry optimization and 
lead to cr bonded complexes with significant binding energies. 
These examples include Figs. 4-6. Open shell dimer cations 
of ethylene and H 2 CNH 2 , and the complex of ethylene with 
H 2 CNH + make stable and strong n complexes, Figs. 7 and 8. 

To summarize, three types of complexes are predicted to 
form from the simple association of ions and molecules included 
in this study. Strong covalently bound molecules, ion-dipole 
bound complexes, and novel ^-complexes of open shell cations 
are predicted. Most of the cationic species considered here have 
large oscillator strength electronic transitions in both the near 
and vacuum ultraviolet, and sometimes even in the visible region 
of the electromagnetic spectrum. Many of the electronic transi- 
tions with large oscillator strengths are due to charge- transfer 
type excitations between the monomer units, and are thus likely 
to be easily observed. The open- shell ^-complexes, interestingly, 
have very large oscillator strengths in the visible region of the 
electromagnetic spectrum. This makes them possible candidates 
for some of the interesting observed interstellar emission fea- 
tures that extend throughout the UV-visible and well into the 
infrared electromagnetic spectrum. 

All of the complexes studied here are formed by ion- 
molecule reactions without an association barrier. This point 
is supported by test ion-molecular trajectory calculations 
acetylene* acetylene and ethylene^, which will be reported in 
due course (Bera et al., in prep.). Thus, the question of which 
isomeric structure of a given complex is initially formed depends 
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significantly on the encounter geometry of the ion-molecule 
species, though if there is a large difference in binding energies, 
there can be thermodynamic driving force towards the more sta- 
ble complex. This is a question that we have not addressed in 
this study, but it would be interesting and desirable to perform 
ab initio trajectory calculations for a range of different impact 
factors and relative orientations for the complexes studied here to 
explore the initial isomer distributions. Subsequently, there can 
be internal rearrangements from one isomeric form to another, 
which are controlled by the excess internal energy and potential 
energy surface barriers. Hence, it would be useful to characterize 
the barriers to internal rearrangements in future work. 
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Appendix A: Cartesian co-ordinates (Angstrom) 
of the molecules, ions and complexes - reported 
here at the MP2/aug-cc-pVTZ level of theory 
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Fig. A.l. a) C 2 H 4 -HNCH + . b) C 2 H 4 -HNCH + . c) C 2 H 4 -HNCH\ Boat, 
d) C 2 H 4 -HNCH + , Bridge. 
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Fig.A.2. a) C 2 H 4 

-H 2 CNHJ. b) C 2 H 4 

-h 2 cnhj. 
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Fig. A.3. a) HCCH-HCNH+ . b) HCCH-HNCH+ . c) HCCH-HCHN+, 
Boat. 
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Fig.A.4. a) HCN-C 2 HL 
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7. a) H 2 CCH 2 -H 2 CNH + . 
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5. a) HCCNCH+, Branched, b) HCCNCH\ ] 
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Fig. A.8. a) 

[H 2 CNH 2 ]J. b) [H 2 CNH 2 ]J. 
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Fig. A.6. a) HCCC(H)NC + , Branched, b) HCCCNCH + , Linear. 
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